Orbital currents, anapoles, and magnetic quadrupoles in CuO 



S. Di Matteo 1 and M. R. Norman 2 

1 Groupe theorie, Departement Materiaux et Nanosciences, 
Institut de Physique de Rennes UMR UR1-CNRS 6251, 
Universite de Rennes 1, F-35042 Rennes Cedex, France 
2 Materials Science Division, Argonne National Laboratory, Argonne, IL 60439, USA 

(Dated: March 4, 2013) 

We show that orbital currents in a CuCb plane, if present, should be described by two inde- 
pendent parity and time-reversal odd order parameters, a toroidal dipole (anapole) and a magnetic 
quadrupole. Based on this, we derive the resonant X-ray diffraction cross-section for monoclinic CuO 
at the antiferromagnetic wavevector and show that the two order parameters can be disentangled. 
From our analysis, we examine a recent claim of detecting anapoles in CuO. 

PACS numbers: 78.70.Ck, 75.25.-j, 74.72.Cj 



I. INTRODUCTION 

A theory of the pseudogap phase in high te mpe rature 
superconducting cuprates, proposed by VarmsP^, intro- 
duced intra-unit cell orbital currents within Cu02 planes 
as a key feature characterizing the physical properties 
of this phase. Varma's theory is based on a mean-field 
solution of a three band Hubbard model, and predicted 
two possible orbital current patterns, with their resulting 
orbital moments pointing out of the plane. Photoemis- 
sion dichroism 3 and neutron scattering experiments^"^ 
are consistent with one of these patterns, depicted in 
Fig. [T] This current pattern can be characterized by an 
in-plane vector order parameter, the toroidal moment (or 
anapole)P. This vector, odd under inversion and time- 
reversal, has the same symmetry properties as the mag- 
netoelectric susceptibility, which it is usuall^ identified 
with. For the Cu02 case, the anapole has two possible in- 
plane orientations 9 , forming an E n representation. But 
the neutron experiments 4 in addition show that the ob- 
served magnetic moments have an in-plane component, 
which is not consistent with solely in-plane currents if 
they have an orbital origin. Very recently, a resonant X- 
ray diffraction (RXD) measurement performed at the Cu 
L3 edge also claimed the detection of orbital currents in 
the commensurate antiferromagnetic phase of monoclinic 
CuO, where the inferred anapole vector does not lie in 
the Cu02 planes^. 

The aims of the present article are the following: (a) 
for a proper characterization of the orbital currents, a 
second order parameter, the magnetic quadrupole, is 
needed; (b) the direction of the inferred toroidal moment 
in CuO is actually 90° rotated from that illustrated in 
Ref. [lOj (c) a toroidal origin for the observed RXD signal 
at the Cu L3 edge in CuO is unlikely since the El-Ml 
matrix elements are extremely small - rather, such ef- 
fects would be more visible at the Cu K edge through 
E1-E2 interference. And although symmetry allowed, we 
believe anapoles are more likely in the higher temper- 
ature multiferroic 11 phase of CuO, where the moments 
form a helhP^. Regardless, we emphasize the presence 
of a second, independent, order parameter, the magnetic 



quadrupole, that is time-reversal and parity odd, and has 
in principle the same order of magnitude as the toroidal 
dipole, and therefore should play a role in any description 
of orbital currents. Though the presence of the magnetic 
quadrupole does not change the theory of Ref. [9] qualita- 
tively, some details have to be reconsidered. We conclude 
by highlighting some general ambiguities in the detection 
of the orbital currents by means of resonant X-ray diffrac- 
tion. 



II. TOROIDAL DIPOLE AND MAGNETIC 
QUADRUPOLE 

The toroidal dipole and magnetic quadrupole are inti- 
mately connected, being, respectively, the antisymmetric 
and the traceless symmetric parts of the bilinear form 
r a mp, where r is the position vector of the magnetic 
moment fa = /ib(£ + 2S), with \±b = ^ the Bohr 
magneton in SI units and a, f3 = x,y,z. They can 
be defined from the expansion of the magnetic energy 
Wm = J dfJ(r) • A(f) up to second order in the field 
(here J(r) is the current density and A(f) t he vect or po- 
tential at position r). This can be written as^HI] W m = 
-m-B+\Y;ij Mij fi(0)+f-VxB. Besides the mag- 
netic moment m, we defined the magnetic quadrupole 
= \ j dV [3(xirrij + rriiXj) — 2{x • fn)5ij] and the 
magnetic (polar) toroidal moment: T = J dV(x x mj^. 

From this, it is possible to identify the toroidal dipole 
T as the second order term (equivalent order as the mag- 
netic quadrupole) that couples to the curl of the magnetic 
field (i.e., to currents and displacement currents). It is 
worthwhile to underline this result in the light of recent 
literatur e! 8 1 141 : the presence of toroidal dipoles T does not 
necessarily imply a magnetoelectric coupling of the kind 
T • (E x B). A symmetry analysis alone does not allow 
one to decide whether the T-dipole is rather coupled to 
Vx5, which is clearly not magnetoelectric. This means 
that in a magnetoelectric material with an antisymmetric 
magnetoelectric tensor otij = —OLj%, one is not allowed to 
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identify the 3 components of this latter with the 3 com- 
ponents of T. In fact, though (V x B) couples with T, we 
are not guaranteed that (E x B) couples with the same 
vector, i.e., with a current loop around a torus. There- 
fore, a material characterized by both a parity odd and 
a time-reversal odd symmetry does not necessarily cou- 
ple to the vector (E x B) in its free-energy expansion. 
A coupling with (V x B) respects these symmetries as 
well, without showing magnetoelectricity. As such, or- 
bital currents, that can be described through T and M^ , 
are not necessarily associated with magnetoelectricity. 

In general, whatever configuration of two oppositely 
oriented orbital currents occurs, it is always characterized 
by a magnetic quadrupole that can be associated with 
their flow, as pictorially shown in Figs, [I] and [2] 

Consider a reference coordinate system xyz and a mag- 
netic distribution whose total net magnetization is zero, 
like the ones depicted in Figs. [I] and [5] (two magnetic mo- 
ments m a at positions r**, a = 1, 2). It is possible to de- 
fine the magnetic quadrupole and the magnetic toroidal 
dipole for these configurations as the first magnetic mul- 
tipoles that are different from zero 
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The magnetic quadrupole is the traceless symmetric 
part of the cartesian tensor m^- and the toroidal dipole 
is the antisymmetric part (we can then call T z = T xy , 
T y = T zx and T x = T yz ). 

We can evaluate and for the orbital currents of 
Fig.[l] The two loop currents are equivalent, by Ampere's 
theorem, to two magnetic moments m directed along ±z 
at positions (—1,1,0) and (1,-1,0), respectively. We 
get therefore M xz = —M yz — T x = T y . It should be 
noticed that by rotating the axes of Fig. [1] by 45°, only 
one component of M and one of T differ from zero. This 
result is independent of the choice of the origin, as the 
toroidal dipole and the magnetic quadrupole are the first 
non-zero multipoles of this magnetic distribution. 

One might ask whether, by continuously deforming the 
relative positions of the two current loops, it is possible 
to find a configuration where the magnetic quadrupole is 
zero and only the toroidal dipole is present or vice versa. 
The configuration of Fig. [I] is the only one having an 
equal value of toroidal dipole and magnetic quadrupole 
components, and any rotation of the magnetic dipoles 
(with the constraint that the total magnetic moment of 
the configuration is zero) leads to an increase of and 
a decrease of T. In the extreme case represented by the 
configuration of Fig. [2j the toroidal dipole is zero. No- 
tice that there is no way to get a pure anapole from a 
two loop current. Indeed, it is not possible to get a pure 
anapole from any n-loop current, unless n goes to infin- 
ity and one gets a pure circulation of magnetic moments 
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FIG. 1: (Color online) Two loop current pattern for Cu02 
planes 1 -^. The magnetic moment fa at f— (-1,1,0) is in the -\-z 
direction and the magnetic moment at -r=( 1,-1,0) is in the 
-z direction. This configuration, characterized by the bilinear 
Txvn z where x is along (-1,1,0), is an equal admixture of a 
toroidal dipole and a magnetic quadrupole. 



FIG. 2: (Color online) Two loop current pattern for Cu02 
planes: the pure magnetic quadrupole case. The magnetic 
moment m at f=(- 1,1,0) is in the xy-plane with components 
(— ra/v^, ?ti/a/2, 0) and the magnetic moment at - f—{ 1,-1,0) 
is in the x^/-plane with components (ra/v 7 ^, — ra/v 7 ^, 0). This 
configuration, characterized by the bilinear rjTOj where x is 
along (-1,1,0), is a pure magnetic quadrupole. 
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(equivalent to a current flowing around a torus). These 
considerations are important in order to build the cor- 
rect order parameter for the two-loop current pattern of 
Cu02 planes, by allowing out-of-plane components of the 
currents. In fact, the relative weight of and is a 
measure of the relative orientation of the current loops, 
as shown in two extreme cases in Figs. [1] and [2] One way 
to measure both order parameters is by means of RXD. 

In the next section, we apply these considerations to 
the case of monoclinic CuO. 



III. CRITICAL ANALYSIS OF RXD IN CUO 

Scagnoli et al. measured the Q=(l/2, 0,-1/2) magnetic 
Bragg reflection in the low temperature monoclinic (Cc) 
phase of CuO at the Cu L3 edge (i.e., in the commensu- 
rate antiferromagnetic phase). They found a dichroic sig- 
nal with a sin(2^) dependence, where ip is the azimuthal 
angle about Q. They interpreted their signal as evidence 
of El-Ml (electric dipole-magnetic dipole) interference, 
which they related to the presence of a toroidal moment. 
We shall analyze this claim by repeating the calculations 
using a somewhat more transparent formalism based on 
cartesian tensors. 

We begin with a description of the CuO commen- 




surate antiferromagnetic structure, which is shown in 
Fig. 3. This structure consists of (1,0,-1) antiferromag- 
netic chains which are intersected by (1,0,1) ferromag- 
netic chains. Ribbons of Cu02 planes run along the 
(1,1,0) and (-1,1,0) directions. We describe the mag- 
netic unit cell using 16 copper ions with doubled a and 
c lattice constants; therefore the (1/2,0, —1/2) reflection 
becomes (1,0,-1) in our notation. This unit cell corre- 
sponds to the four blocks shown in Fig. 4. Copper ions 
are at the (equivalent) positions: Cui = (1/8,2/, 0) - (E) 
Cu 2 = (3/8,2/ + 1/2,0) - (£); Cu 3 = (5/8, y, 0) - (f) 
Cu 4 = (7/8,j/+l/2,0) - (T); Cu 5 = (1/8, y, 1/4) - (fm b ) 
Cu 6 = (3/8,S7+l/2, 1/4) - (m 6 ); Cu 7 = (5/8,17,1/4) - 
(m 6 ); Cu 8 = (7/8,17+1/2,1/4) - (fm h ). In parenthe- 
sis, we included the symmetry with respect to the first 
ion Cui: E is the identity, T is time-reversal, 7715 is the 
glide plane perpendicular to the 6-axis, and Tm^ is the 
product of the last two operations. The 8 other ions are 
obtained by a translation of (0,0, 1/2), multiplied by T, 
as is clear from Fig. 4. 

The structure factor is defined as F(hkl) = 
J2j=i fj{uj)e 2 ' K ' L ^ hx ^ Jrkv ^ +lz 3\ where Xj, yj and Zj are the 
three fractional coordinates for each of the 16 ions in the 
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FIG. 3: (Color online) Commensurate antiferromagnetic 
structure of monoclinic CuO. Small (blue) atoms are cop- 
per, large (red) oxygen, with arrows denoting the magnetic 
moments pointing along ±6. The antiferromagnetic chains 
run along (1,0,-1) (lower left to upper right; green), with in- 
tersecting (1,0,1) ferromagnetic chains (upper left to lower 
right; blue). Ribbons of CuCb planes run along the (1,1,0) 
and (-1,1,0) directions. 



FIG. 4: (Color online) Commensurate antiferromagnetic CuO 
unit cell, equal to four times the crystallographic cell, pro- 
jected in the ac-plane. The fractional b coordinate of each ion 
is explicitly written. Open (blue) circles represent Cu ions 
having a magnetic moment pointing along b, closed (blue) cir- 
cles along -6, with oxygen ions in the alternate rows (closed 
red circles). The anapoles (arrows) lie in the ac plane, or- 
thogonal to Q. 
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unit cell denned above, and fj(u) is the dynamical atomic 
scattering amplitude (u is the photon energy) defined in 
Eq. |4j Here (hkl) = (1,0,-1). As all 16 Cu ions are 
related by some symmetry elements, the structure fac- 
tor can be reduced to the dependence of only one atomic 
scattering factoi^E^, say that of j\ (Cui). We get 



2e^/ 4 (l + i)(l-T)(l + m 6 )A 



(3) 



The final result can be read as follows: the total am- 
plitude (which must be squared to get the intensity) is 
provided by the time-reversal odd (magnetic, because of 
the term (1 — T)) part of the atomic scattering factor 
of Cui that is also even under the glide plane (were 
it odd, the last factor (1 + m&) would have given zero). 
This result is important because it rules out a helical 
magnetic pattern arising from any residual component of 
the higher temperature multiferroic phase 18 as a cause 
of the signal. In fact, only the component of the mag- 
netic moment along the b axis is detectable at the (1,0,-1) 
reflection, as the a and c components are odd under m&. 

The atomic scattering factor f\ appearing in Eq. |3J 
neglecting for the moment E2 (electric quadrupole) pro- 
cesses and considering as in Ref. 10 just El and Ml terms, 
can be written as h = A E1 " E1 + / E1_M1 + /i M1 " E1 + 
yMi-Mi rpj^ orc [ er Q f magnitude of these terms has 
been calculated by means of the FDMNES program^ 
and Ml matrix elements have been found to be less 
than four orders of magnitude smaller than the El ma- 
trix elements. Therefore we can definitively neglect the 
jMi-Mi am plitudes and the measured signal is therefore 

T „ I /El -El | /El-Ml | /Ml-Eli2 

Ik \ fi +/i +/i I - 
We define the above quantities as 



/iE i-ei = £ A n <* fl |6; • r1* n )<* n |6-; • fl* fl > (4) 



where the resonant denominator A n = {hw — (E n — 
Eg)) -1 provides the uj dependence (here we neglected 
damping, not important for symmetry considerations) 
and e ^ are the outgoing and incoming polarizations, 
respectively. It is common practice to expand the two 
scalar products and then factorize the terms depending 
on the electromagnetic wave (the polarization in the El- 
El case), and those depending on matter. This is usually 
done in spherical tensors^Hl, but it can be done in carte- 
sian tensors as welffl. We get 



A e1 - e1 = £t q/3 f ( 
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where T a p = {{*)&{£%) p represents the electromagnetic 
wave and F ai3 = J2 n A n (Vg\r a \V n )($ n \rp\V g ) repre- 
sents the properties of the sample, with a, j3 = x, y, z. 
With the same notation, we can write 



/r- Mi +/r- Ei =E A « ( 6 ) 

n 

i^gK ■ r\^n)(^n\(ki x el) • (L + 2S)\* g ) 
+ (V g \{k x e* Q ) ■ (L + 2S)\y n )(y n \ei-f\y g ) = 
= J2 [{e*o)a(ki x tijpSctP + {ko x e* o ) a (ei) S afz 

a,p 

= Yl [(£*o)a(ki X Si) p (tftS a p + i%S a p) 

a, 13 

+ (ko x Z) a (ei)f)(3tSp a - i%Sp a ) 
= [MS a p((? ) a (ki x ei) e + (ei) a (k x 



where we denned the two hermitian conjugate quanti- 
ties S al3 = E„A„(* g |(r)a|*n)(*n|(i + 25) /3 |* g ) and 

Sap = En A n(^ g \{L + 2S) a \^ n )^ n \(r)^ g ). Notice 
that for the El-El term, the antisymmetric and complex 
conjugate coincide for F a p, but this is not the same for 

S a p and S a p: ^[S a p] = ^[Sp a ] and $t[S a p] = -$f[Sp a ]. 
This relation was used in the next to last step of Eq. [6] 
There are two independent, real 3x3 matrices appear- 
ing in Eq.[6j $tS a ,p, time-reversal odd, and i$sS a ^p, time- 
reversal even, as can be seen from their definitions 



®S° P = ~ E A " [(* fl l(* ! )a|*n><*n|(i + 2^) |* fl 



-(9 g \(L + 2§) J*„><*„|(»Da|* fl > 



(7) 



i$!S a p = -^A n \(* g \(r) a \* n )(* n \(L + 2S) |* fl ) 



-<* fl |(L + 2^)j* n )<* n |(f)a|*fl> 



(8) 



We remind that the time-reversal operation performs a 
complex conjugation and reverses the sign of L and S. 
Notice the difference with the El-El 3x3 tensor, char- 
acterized by just 9 independent components, such that 
the symmetric (real) elements are time-reversal even and 
the antisymmetric (imaginary) elements are time-reversal 
odd. Instead, in the El-Ml case, there are 18 terms, 9 of 
which are time-reversal even (imaginary, independently 
of spatial symmetry) and 9 are time-reversal odd (real, 
again independently of spatial symmetry). With spher- 
ical tensors, one can identify each of the 9 components 
in terms of irreducible tensors as C = (1 component), 
C = 1 (3 components) and C = 2 (5 components). The 
correspondence of these irreducible tensors with carte- 
sian tensors is very simple: C = corresponds to the 
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trace of the matrix, C = 1 is the antisymmetric part, 
C = 2 is the traceless symmetric part. This is valid for 
all the El-El amplitudes and each of the time-reversal 
odd/even matrices of El-Ml events. 

In the light of the experiment on CuO, we focus on 
time-reversal odd terms, $tS a p, and use the following 
spherical linear combinations: 



[$ts} {0) = J2^s Q 



[$tS] 



(i) 



2 [® S aP 



[fits] 



(2) _ 1 



ytS a /3 + $tSp a — -SaffiSap 



(9) 



(10) 
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The pseudoscalar part (Eq. [9| is not effective in our 
case as it is glide-plane odd. The antisymmetric term 
(Eq. [To]) is the anapole (polar toroidal dipole), whereas 
the symmetric, rank- two tensor (Eq. [TT| ) is the magnetic 
quadrupole. 

For actual calculations, we choose a reference frame 
with the z direction along Q = (1, 0, —1), the y direction 
along the b axis, and the x direction along the (—1, 0, —1) 
direction (real space), in order to have a right-handed 
frame. In our frame, the sample is fixed and the x-ray 
beam rotates clockwise, in order to match the experi- 
mental condition where the beam is fixed and the sam- 
ple rotates counterclockwise. According to the choice 
of Scagnoli et al, the azimuthal angle ip is zero when 
the diffraction plane contains the b monoclinic axis (see 
Fig. [5| . An intrinsic ambiguity of sign (both ip and ip + 7r 
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FIG. 5: Diffraction plane defined by the vectors and k > It 
is drawn for -0 = 0, when the diffraction plane contains the b 
monoclinic axis. 



satisfy this condition) does not affect the final results 
for the anapole, but can affect the sign of the magnetic 
quadrupole contribution (see Eq. 18 below). Some defini- 
tions useful for the following, in our xyz frame of Fig. [5j 
are: 

4j e oa = (-cos^,sin^,0) 
e i7T = (sin sin ip : sin cos ip, cos 0) 
ton — (— sin # sin?/;, — sin # cos?/;, cos#) 
ki = k(cos0s'mip, cos # cos?/;, — sin (9) 
k Q = /c(cos#sin?/;, cos#cos?/;, sin#) 
Q k Q -ki = 2/c(0,0,sin<9) 
ki + k Q = 2/c(cos#sin?/;, cos # cos?/;, 0) 

where k is the modulus of the wavevector and is the 
Bragg angle. 

The El-El atomic scattering factor is known to be^ 



fi 



El-El 



■ m 



(12) 



Instead, the El-Ml atomic scattering amplitudes con- 
tributing to the structure factor (3j) are the components 
of the anapole in the xz plane and the xy and yz compo- 
nents of the magnetic quadrupole. The reason is that the 
anapole is a time-reversal odd, polar vector, and there- 
fore its y component is odd under m^, contrary to the 
x and z components. The magnetic quadrupole is such 
that only the xz and yz components are non zero after 
the application of 1 + rhb in Eq. [3j The contribution pro- 
portional to the anapole is the vector term from Eq. [lOj 
whose polarization and wavevector components are given 
by P = (e£) x (ki x e-) - (k Q x e£) x (e-). Therefore Eq.g 

limited to this term, becomes: fa cx T • P, where T 
is the toroidal dipole (anapole) vector^. For the mag- 
netic quadrupole, similar rules apply. In particular, for 
the M rs component, we have fa cx M rs (^t ) r (ki x e^) s + 
(€* ) s (ki x €i) r + (ei) r (k x e^) s + (e±) 3 (k x e£) r ). 

From these definitions we can evaluate the scattering 
amplitudes in the a, 7r and circular ± polarization con- 
ditions. The El-El magnetic structure factors are 

f™-*\aa) = 

/f 1_E1 ((T7r) cx —irriy cos cosip 



ff 1 E1 (7TCr) CX IVfly COS COS Ip 

/f 1_E1 (7r7r) cx im y sm(20)smip 

whereas those corresponding to fa are: 

fa( <J(J ) (X- %T X cos sin ip 

fa(cnr) cx -T x sin(20) cos ip 

/ia(7rcr) cx T x sm(20) cos ip 

fai^Tr) cx 2T X cos sin ip 



(13) 



(14) 
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and those corresponding to fib are: 

fibber) oc 2M yz cos 9 sin ip 

fib(°"K) oc M xy cos 2 6 sin(2 , 0) 

fibber) oc M xy cos 2 sm(2ip) 

fib(7T7r) oc — 2M yz cos sin ip 



(15) 



Amplitudes for circular polarizations are defined as 
= f aa ± if a7T . Circular dichroism can be defined 
in two ways, according to whether one measures or not 
the outgoing polarizations. In the former case if (say) 
outgoing measured polarization is a, one should evaluate 
AI a = 1+^ — I- y(J . If no outgoing polarization is mea- 
sured, then AI = + J +?7r — I _ a — J_ n should be 
evaluated. The latter case is that of Ref. From the 

previous definitions we get: 



AI = |F + .J 2 + |F +)OT | 2 -|F_, CT | 2 -|F_ )7r | 2 (16) 



^ \f(J,(T H~ ^flT,a\ H~ \fa,7T H~ ^«/*7T,7T | 
|/cr,CT ^flT,a\ |«/ct,7T ^ «/*7T , 7T | 

CX- ^-fa,a^fiv,a ^ fa ,iv ^ fir, it 
H~ ^.fir,a^fa,a ^ fir , it ^ fa, it 

Analogously: 

AI CT = |F + , a | 2 - |F_, a | 2 oc + &U,*$Sf*,* 

(17) 

Taking the imaginary quantities from Eq. [l3j and the 



real quantities from Eqs. 14 and 15, we get for the first 
dichroism: 



AI oc m y T x cos 4 sin(2^) 

+ m y M xy cos 3 sin sin(2 , 0) sin ip (18) 

and for the second: 



AI a oc m y (T x + M^) cos 2 6> sin(2^) 



(19) 



Therefore from a coupled analysis of the azimuthal 
scan of both kinds of dichroism, it is possible to obtain 
the relative value of all three components T x , M xy and 
M yz : the simp correction to the sin(2^) azimuthal scan 
is proportional to M xy (if the values of m y and of the 
radial transition matrix elements are known). More in- 
teresting is the information that we can get from the ratio 
of the sin(2?/0 modulation of AI and AI a : 4AI a /AI = 
cos -2 0(1 + M yz /T x ), which is an absolute measurement 
of the magnetic quadrupole to toroidal dipole ratio, in- 
dependent of m y and of any radial transition matrix el- 
ement. This would be important in the light of the find- 
ings of Section II, concerning the ratio of the anapole 
to magnetic quadrupole as a measure of the out-of-plane 
component of the orbital currents. 

We also performed numerical calculations by means of 
the FDMNES program at the Cu pre K-edge at the same 



reflection. The aim is to find out whether clearer evidence 
of orbital currents can be obtained through E1-E2 ma- 
trix elements, e.g., transitions from Is to Ap levels and 
then from 3d back to Is, by taking advantage of the pd 
hybridization in CuO. The calculated order of magnitude 
of the signal in the era channel, though small ~ 4-10 _4 rQ, 
is well within the present sensitivity at synchrotron facili- 
ties. However, several multipole components are present, 
as detailed below, and in order to identify each term, 
one should take advantage of both azimuthal scans and 
measurements of several reflections of the same class as 
(1, 0, —1), in order to change the value of as well. The 
measured components at the Cu pre K-edge are the fol- 
lowing: in the E1-E2 channel we get T x , T z , M xy , M yz , 
7>, T xz 2, T xyz , T x ( x 2_ 3y 2^ where T a ^ is the toroidal 
octupole. Moreover, at the pre K-edge, E2-E2 terms in 
the a a channel are present, of the same order of mag- 
nitude (~ 10 _4 ro). They are: m y: O yz 2, O z ( x 2_ y 2^ and 
O y (3 X 2- y 2), where O a f3 1 is the magnetic octupole. Focus- 
ing on just the toroidal dipole and magnetic quadrupole, 
their azimuthal dependence is the followin g 25 * 26 !: 



fx (^cr) ex 2 cos 6 sin?/; 



/x(<j7r) (x sin(2#)cos^ 



T x - 2 -^M yz 



cos 2 0sm(2ip)M xy 



fie (kg) (x — sin (20) cosip 



T x - ^M yz 



cos 2 <9sin(2^)M^ 



fx (tttt) oc 2 cos 6 sin ip 

(cos 2 (9-3 sin 2 0) 



T x - ^yM yz 



(20) 



This could be compared to future experiments. In- 
terestingly, the full structure factor for a generic (h, 0, 1) 
reflection reads: 



F(hOl) = e™/\l + (-) h f)(l + (-) l f)(l + (-i) h f) 
(1 + (i) h+l m h )h (21) 

an expression that reduces to Eq. ([3| when h = 1 and 
I = — 1. This shows that magnetic reflections are charac- 
terized by both h odd and I odd. The absence of mag- 
netic reflection as measured by neutron scattering for 
h + I = An + 2 showed that the magnetic moment has 
only the y component (otherwise the glide plane would 
have implied a canted x or z component that would have 
given an intensity to this class of reflections). However, 
this extinction rule, implied by the (1 + (i)^ 1 ^) term, 
might be violated for E1-E2 terms. A signal at, e.g., 
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the (1, 0, 1) reflection is in this case a direct indication of 
a magnetic quadrupole or toroidal dipole and should be 
looked for. 



IV. CONCLUSIONS 

Given the smallness of the Ml matrix elements, we find 
that the dichroism signal observed in Ref.fTQlis unlikely to 
be due to El-Ml interference. And although anapoles are 
symmetry allowed for the commensurate antiferromag- 
netic phase of CuO, we find them physically unlikely, as 
the displacement of the Cu ions off of their C2/c locations 
in the Cc structure is along 6^3 which is the direction of 
the magnetic moments. The fact that the resulting cross 
product is zero casts doubts on the existence of anapoles 
arising from a magnetoelectric effect. 

Recently, Joly has suggested that the aa signal 
observed by Scagnoli et air** is most likely due to 
birefringence^. That is, because of the biaxial nature 
of the monoclinic structure, the light will rotate in the 
sample, causing a an signal of purely magnetic dipole ori- 
gin to show up in the a a channel. In fact, the split nature 



of the RXD peak observed in Ref. [10] is almost certainly 
due to strong self-absorption, leading further credence 
to this scenario. Whether such a scenario can account 
for the pronounced sin(2^) azimuthal dependence of the 
observed dichroism signal remains to be seen. 

From the above considerations, it would be extremely 
interesting to look at RXD in the multiferroic phase, es- 
pecially at the pre K-edge, where birefringence effects 
should be much less important. In this phase, the mag- 
netic moments rotate in a plane that contains b and a vec- 
tor near (1,0, 3) 12 , with a ferroelectric polarization along 
Zffl Moreover, such experiments in the cuprates would 
be very illuminating in resolving the question of whether 
an orbital current description is appropriate for the ob- 
served neutron scattering signal in the pseudogap phase. 
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